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Abstract 

The instability against emission of massless particles by the trapping horizon of an evolving black 
hole is analyzed with the use of the Hamilton- Jacobi method. The method automatically selects one 



special expression for the surface gravity of a changing horizon. Indeed, the strength of the horizon 



singularity turns out to be governed by the surface gravity as was defined a decade ago by Hayward 
using Kodama's theory of spherically symmetric gravitational fields. The theory also applies to point 
masses embedded in an expanding universe, were the surface gravity is still related to Kodama- 
Hayward theory. As a bonus of the tunneling method, we gain the insight that the surface gravity 
still defines a temperature parameter as long as the evolution is sufficiently slow that the black hole 

Oh' 

■ pass through a sequence of quasi-equilibrium states. 
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semi-classical theory of black hole radiation and evaporation has perhaps reached a satisfactory state of 



With the exception of few lower dimensional models where dynamical computations can be done, the 



development only for the case of stationary black holes. It includes now a vast range of topics, ranging from 
the various derivations of the quantum Hawking's effect to the efforts aiming to a statistical interpretation 



of the area law and the associated thermodynamical description. One of the surprising aspects of these 
results is that the radiation caused by the changing metric of the collapsing star approaches a steady 
outgoing flow for large times, implying a drastic violation of energy conservation if one neglects the 
back reaction of the quantum radiation on the causal structure of spacetime. But the back reaction 
problem has not been solved yet in a satisfactory way. As pointed out by Fredenhagen and Haag long 
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ago [1], if the back reaction is taken into account by letting the mass of the black hole to change with 
time, then the radiation will originate from the surface of the black hole at all times after its formation 
and "it will no longer be precisely calculable from a scaling limit" [1] (the term originates from the 
link between the Hawking's radiation and the short distance behaviour on the horizon of the two-point 
function of the quantum field). Thus one is concerned to show, in the first place, that at least some 
sort of instability really occurs near the horizon of the changing black hole. This question is non trivial 
since a changing horizon is not even a null hypersurface, although it is still one of infinite red shift. In 
this vein we shall analyze this question for a class of dynamical black hole solutions that was inspired 
by problems not directly related to black hole physics, although these were subsequently reconsidered in 
the light of the black hole back reaction problem in the early Eighties. The metrics we shall consider 
are the Vaidya radiating metric [2], as revisited by J. Bardeen [3] and J. York [4], together with what 
really is a fake dynamical black holes, the McVittie solution representing, in author's mind, a point mass 
in cosmology [5]. And the strategy we shall use is a variant of a by now well known method due to 
F. Wilczck and M. Parikh [6] according to which the process of Hawking radiation is akin to a tunneling 
effect through the horizon. The method was refined and extended to more general cases in [7] and others 
papers as well [8]. For criticism and counter criticism see also [9]. It must be kept in mind that these 
solutions do not correspond precisely to the standard notion of a black hole, but they do have horizons 
(apparent and/or trapping) to which the familiar black hole theorems seems to apply [10]. 

The variant referred to above is the Hamilton- Jacobi method introduced in [11], so called after the 
comparison analysis with the Parikh- Wilczek method done by B. Kerner and R. Mann [12]. The tunneling 
method provides not only new physical insight to an understanding of the black hole radiation, but is also a 
powerful way to compute the surface gravity for a vast range of solutions. As a matter of facts, in the past 
decade some different definitions of the surface gravity of an evolving black hole were proposed [10], [13], 
until one which met with all requirements was introduced in [14] using a key ingredient invented by 
H. Kodama [15] (more on this later). The HJ method can also be applied to the more elaborate theory 
of isolated horizons of Ashtekhar and co-workers (see [16] and references therein. The literature is 
quite extensive). Results in this direction appeared recently in [17], where the Parikh- Wilczek and HJ 
methods are compared and showed to agree. An early study of the evolution of evaporating black holes 
in inflationary cosmology is in [18]. 

We begin in Sec. II to present the relevant solutions to which the method will be applied. This will 
actually be done in Sec. Ill, where we shall also display the "generalized first law" and the conclusions 
to be drawn from it. 
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2 Dynamical black holes, and a fake one 

We consider first spherically symmetric spacetimes which outside the horizon (if there is one) are described 
by a metric of the form 



ds 2 = -e™^A(r, v)dv 2 + 2e* {r ' v Uvdr + r 2 dS 2 D „ 



(2.1) 



where the coordinate r is the areal radius commonly used in relation to spherical symmetry and v is 
intended to be an advanced null coordinate. In an asymptotically flat context one can always write (we 
use geometrized units in which the Newton constant G = 1) A(r, v) = 1 — 2m(r, v)/r D ~ 3 . This metric was 
first proposed by Vaidya [2] in D = 4 dimensions, and studied in an interesting paper during the classical 
era of black hole physics by Lindquist et al [19]. It has been generalized to Einstein-Maxwell systems and 
de Sitter space by Bonnor- Vaidya and Mallet, respectively [20]. In the special case <&(r, v) = \P(r, v) it 
was then extensively used by Bardeen [3] and York [4] in their semi-classical analysis of the back reaction 
problem, the former to establish the stability of the event horizon in the geometry modified by Hawking 
radiation back reaction (contrary to a claim of F. Tipler), the latter in an attempt to explain dynamically 
the origin of the entropy of the black holes. We shall call it the Vaidya-Bardeen metric. If one wishes the 
metric can also be written in double-null form. In the (v,r)-plane one can introduce null coordinates x^ 
such that the dynamical Vaidya-Bardeen BHs may be written as 



ds 2 = -2f(x + ,x-)dx + dx~ + r 2 (x+, x~)dS 2 D _ 2 , 



(2.2) 



for some diffcrentiable function /. The remaining angular coordinates contained in dSj)_ 2 do not play 
any essential role. In the following we shall use both forms of the metric, depending on computational 
convenience. The field equations of the Vaidya-Bardeen metric in D = 4 dimensions are of interest. They 
read 

§= = -4"X, (2.3) 

The second example we are interested in is the McVittie solution [5] for a point mass in a Friedmann- 
Robertson- Walker flat cosmology. In D-dimensional spacetime in isotropic spatial coordinates it is given 
by [21] 

ds 2 = -A(p, t)dt 2 + B(p, t) (dp 2 + p 2 dS 2 D _ 2 ) (2.4) 



with 



A(p,t) 



1 - 



a(t)p 



D-3 



1 + 



\a(t)p) 



D-3 



B(p,t) = a(t) 2 



a(t)p 



D-3 



2/(13-3) 



(2.5) 



When the mass parameter m = 0, it reduces to a spatially flat FRW solution with scale factor a(t); when 
a(t) = 1 it reduces to the Schwarzschild metric with mass m. In four dimensions this solution has had a 
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strong impact on the general problem of matching the Schwarzschild solution with cosmology, a problem 
faced also by Einstein and Dirac. Besides McVittie, it has been extensively studied by Nolan in a series 
of papers [22] . To put the metric in the general form of Kodama theory, we use what may be called the 
Nolan gauge, in which the metric reads 

ds 2 =-(A a - H 2 {t)r 2 ) dt 2 + A^dr 2 - 2A- 1/2 H(t)r drdt + r 2 dS 2 D _ 2 (2.6) 

where H(t) = dja is the Hubble parameter and, for example, in the charged 4-dimensional case, A s — 
1 — 2m /r + q 2 /r 2 and in D dimension A s = 1 — 2Al/r D ~ 3 + Q 2 /r 2D ~ 6 . In passing to the Nolan gauge 
a choice of sign in the cross term drdt has been done, corresponding to an expanding universe; the 
transformation H(t) — » —H(t) changes this into a contracting one. In the following we shall consider 
D = 4 and q = 0; then the Einstcin-Friedmann equations read 

3H 2 = 8np , 2A- 1/2 H{t) + 3H 2 = -87175 . (2.7) 

It follows that A s — 0, or r = 2m, is a curvature singularity. In fact, it plays the role that r = has 
in FRW models, namely it is a big bang singularity. When H = one has the Schwarzschild solution. 
Note how the term H 2 r 2 in the metric strongly resembles a varying cosmological constant; in fact for 
H a constant, it reduces to the Schwarzschild-de Sitter solution in Painlcvc coordinates. As we will see, 
the McVittie solution possesses in general both apparent and trapping horizons, and the spacetime is 
dynamical. However, it is really not a dynamical black hole in the sense we used it above, since the 
mass parameter is strictly constant: for this reason we called it a fake dynamical BH. This observation 
prompts one immediately for an obvious extension of the solution: to replace the mass parameter by a 
function of time and radius, but this will not be pursued here. 

The solutions being given, one may ask whether they have anything to do with black holes. Unfor- 
tunately it is not entirely clear what should be considered a black hole in a dynamical regime. However, 
some kind of horizon must be present, and moreover, the evolution should be sufficiently slow to per- 
mit a comparison with the more familiar stationary case. In order for a horizon to exist some metric 
components must have a zero somewhere. To clarify the issue we have to introduce few definitions, so 
at this point we have to refer the reader to references [10], [14], [23], [24] and [25] for more details. The 
expansions 8± of the null geodesic congruences orthogonal to a closed surface S with codimension two 
and measure fi(S) = flrj-2r D ~ 2 , are defined by 6± = fi(S)~ 1 d±n(S) = (D — 2)r~ 1 d±r. A marginal 
surface S is a (D — 2)-dimcnsional spacelike surface with vanishing expansion, and a dynamical horizon is 
a hypersurface which is foliated by marginal surfaces. It is possible to show that the area of such marginal 
surfaces is non decreasing if the energy inflow is non negative. A future dynamical horizon, say H + , is 
the hypersurface implicitly determined by the condition 9 + = It will be called a trapping horizon if 



4 There are dual definitions involving 0„ and some interchanging of "future" with "past" , see the references for details. 
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+ is strictly decreasing on crossing the horizon from the outside (so that d-6 + < on H + ), a condition 
which insures the non vanishing of the surface gravity (indeed it is violated for extremal black holes and 
naked singularities). We just mentioned the surface gravity; a geometrical definition of this quantity for 
a trapping horizon is [14], n = (D — 2)~ 2 g H \(D — 2)d-6 + + 6 + 6-]^ e Q . We stress that this quantity 
is not the same surface gravity that was defined few years before this new proposal, although it fitted 
equally well to a generalized first law. Later on we will show that k fixes the expansion of the metric 
near the horizon along a future null direction. 

These definitions look somewhat artificial, but in fact they are very natural and connected directly 
with what is known for the stationary black holes. To see this one notes, following Kodama [15], that 
any metric like (|2.ip or (|2.2p . admits a unique (up to normalization) vector field K a such that K a G a b 
is divergence free, where G a b is the Einstein tensor; for instance, using the double-null form, one finds 
K = — g H (d + rd- — cLr9+), while using (|2.ip one has K — e~®d v . In any case the defining property of K 
show that it is a natural generalization of the time translation Killing field of a static black hole. Moreover, 
by Einstein equations K a T ab will be conserved so for such metrics there exists a natural (Kodama wrote 
"preferable") localizable energy flux and its conservation law. Now consider the expression K a V\bKa,] : 
it is not hard to see that on H + it is proportional to Kb- The proportionality factor, a function in fact, 
is the surface gravity: K a V\bK a -\ = —kK\>. For a Killing vector field VbK a is anti-symmetric so the 
definition reduces to the usual one. The study of black holes requires also a notion of energy; the natural 
choice would be to use the charge associated to Kodama conservation law, but this turns out to be the 
Misncr-Sharp energy, which for a sphere with areal radius r is the same as the Hawking mass [26] (here 
D = 4), E = 2~ 1 r(l — 2~ 1 r 2 g H + Q-). If this definition is found not particularly illuminating, then 
one can use the metric (|2.ip to obtain the equivalent expression g^d^rd v r = 1 — 2E/r, where g^ is 
the reduced metric in the plane normal to the sphere of symmetry (the plane (v,r)). In this form it is 
clearly a generalization of the Schwarzschild mass. As we said, E is just the charge associated to Kodama 
conservation law; as showed by Hayward [23], in vacuo E is also the Schwarzschild energy, at null infinity 
it is the Bondi-Sachs energy and at spatial infinity it reduces to the ADM mass. 

Let us apply this general theory to the two classes of dynamical BH we have considered. Using 
Eq. ([2lL]), we have 9 + = (D - 2)e 2 *"*A(r, v)/2r . The condition 9 + = leads to A(r H ,v) = 0, which 
defines a curve th = ru(v) giving the location of the apparent horizon; writing the solution in the 
Vaidya-Bardeen form, that is with A(r,v) = 1 — 2m(r,v)/r, the Misner-Sharp energy of the black hole 
is E = m(r^f(w),w), and the horizon will be trapping if m (rjy,u) < 1/2, a prime denoting the radial 
derivative. The surface gravity associated with the Vaidya-Bardeen dynamical horizon is 

A'(r,v) m(r H ,v) m'(r H ,v) 1 m'(r H ,v) 

K (V) = 7> \r=r H = 2 = n • ( 2 -8) 

2 rjj r H 2r H r H 

We see the meaning of the trapping condition: it ensures the positivity of the surface gravity. In the case 
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of McVittie BHs, we obtain 9± — ±(_D — 2)(vA j qp Hr)/2rf±, where the functions f± determine null 



One may compute from this the dual 



coordinates x^ such that dx^ = f± y/A s ± Hrj dt ± A s 1 ^ 2 dr 
derivative helds d±. Thus, the future dynamical horizon defined by 8 + = 0, has a radius which is a root 
of the equation y/A~ s — Hrn , which in turn implies A s = H 2 r 2 H . Hence the horizon radius is a function 
of time. The Misner-Sharp mass and the related surface gravity are 



1 »^aj» - 1 r 
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H m H 



E = m + - H(t) r H , «(«) = - A s (r*) - R\ H - — = - R\ K - — . (2.9) 



2i? 2i? 



Note that £ = rjj/2. In the static cases everything agrees with the standard results. The surface gravity 
has an interesting expression in terms of the sources of Einstein equations and the Misner-Sharp mass. Let 
T(id) be the reduced trace of the stress tensor in the space normal to the sphere of symmetry evaluated 
on the horizon H + . For the Vaidya-Bardeen metric it is, by Einstein's equations 

rp _r r v , T r _ 1 9m 
T im-T v +T r --^——^ H 

For the McVittie's solution, this time by Fredmann's equations (|2.7[) one has 

T(2D) = -P + P = 

We have then the formula, k = r^E + 27ITH- T^m. It is worth mentioning the pure FRW case, i.e. 
A s = 1, for which n(t) = — \ H(t) + H/2H^j . We feel that these expressions for the surface gravity are 
non trivial and display deep connections with the emission process. Indeed it is the non vanishing of re 
that is connected with the imaginary part of the action of a massless particle, as we are going to show in 
the next section. 




3 Tunneling within the Hamilton-Jacobi method, and the con- 
clusions 

The essential property of the tunneling method is that the action / of an outgoing massless particle 
emitted from the horizon has an imaginary part which for stationary black holes is Iml = itk,~ 1 E, 
where E is the Killing energy and re the horizon surface gravity. The imaginary part is obtained by 
means of Feynman ie-prescription, as explained in [6, 11]. As a result the particle production rate reads 
r = exp(— 21m/) = exp(— 2ttk~ 1 E) . One then recognizes the Boltzmann factor, from which one deduces 
the well-known temperature Th — k/2tt. But more than this, an explicit expression for re is actually 
obtained in terms of radial derivatives of the metric on the horizon. Let us consider now the case of a 
dynamical black hole in the double-null form. We have for a massless particle along a radial geodesic 
the Hamilton-Jacobi equation d+Id-I = 0. Since the particle is outgoing d—I is not vanishing, and we 
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arrive at the simpler condition d+I = 0. First, let us apply this condition to the Vaidya-Bardeen BH. 
One has then 

2e- xV{r ' v) d v I + A{r, v) d r I = . (3.1) 

Since the particle will move along a future null geodesic, to pick the imaginary part we expand the 
metric along a future null direction starting from an arbitrary event (rff(uo), i>o) on the horizon, i.e. 
A{th{vq), Vq) = 0. Thus, shortening r# (vq) = ro, we have A(r, v) = d r A(ro,vo)Ar + d v A(ro, vq) Av ■ ■ ■ = 
2n(vo)(r — ro) + . . . , since along a null direction at the horizon Av = 0, according to the metric (|2.1|) ; 
here k(vq) is the surface gravity, Eq. (|2.8p . From (|3.1[) and the expansion, d r I has a simple pole at the 
event (tq,vq); as a consequence 

Im / = Im / d r Idr = Ira [ dv , -. 2 £^lM = . (3 . 2 ) 

J J A (ro, vo)(r — tq — i0) k(v ) 

where u)(vq) = e~*^ r °'' Uo ' d v I , is to be identified with the energy of the particle at the time vq. Note that 
the Vaidya-Bardeen metric has a sort of gauge invariance due to conformal reparametrizations of the 
null coordinate v. the map v — > v(v), ty(v,r) — > ^(v,r) + \n(dv/dv) leaves the metric invariant, and the 
energy is gauge invariant too. Thus we see that it is the Hayward-Kodama surface gravity that is relevant 
to the process of particles emission. The emission probability, Y — exp(— 2ttu)(v)/k(v)), has the form of a 
Boltzmann factor, suggesting a locally thermal spectrum. For the McVittie BH, the situation is similar. 



In fact, the condition d+I = becomes d r I = —F(r, t)~ 1 dtl , where F(r, t) = y/ A s (r)(y/ A s (r) — rH(t)). 
As before, we pick the imaginary part by expanding this function at the horizon along a future null 
direction, using the fact that for two neighbouring events on a null direction in the metric (|2.6[) . one has 
t-to = {2H§r^y 1 (r - r ), where H = H(t ). We find the result 

F(r, t) = (l A' s (r ) - r H 2 Q - (r - r ) = «(to)(r - r ) . . . (3.3) 

where this time ro = rjiito). From this equation we see that d r I has a simple pole at the horizon; 
hence, making use again of Feynman ie-prescription, one finds ImJ = 7tk(£o) ^(£o)j where uj(t) = dtl is 
again the energy at time t, in complete agreement with the geometric evaluation of the previous section. 
Obviously, if k vanishes on the horizon there is no simple pole and the black hole should be stable^. 
The kind of instability producing the Hawking flux for stationary black holes evidently persists in the 
dynamical arena, and so long as the evolution is sufficiently slow the black hole seems "to evaporate 
thermally" (paraphrasing [3]). Note that the imaginary part, that is the instability, is attached to the 
horizon all the time, confirming the Frcdcnhagen-Haag suggestion quoted in the introduction. It is worth 
mentioning the role of k in the analogue of the first law for dynamical black holes (contributions to this 
problem for Vaidya black holes were given in [28]). Using the formulas of the projected stress tensor 



3 However, charged extremal black holes can radiate [27]. 
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T(2D) given above, and the expression of the Misner-Sharp energy, one obtains the differential law 

dE=^L-^EldV H (3.4) 

provided all quantities were computed on the horizon. Here Ah = At:r H is the horizon area and Vb = 
Airr^ /3 is a formal horizon volume. If one interprets the "d" operator as a derivative along the future null 
direction one gets Hayward's form of the first law. But one can also interpret the differential operation 
more abstractly, as referring to an ensemble. Indeed, to obtain Eq. Q3.4p it is not necessary to specify 
the meaning of the "d" . It is to be noted that the same law can be proved with other, incquivalcnt 
definitions of the surface gravity, even maintaining the same meaning of the energy. Thus it was not a 
trivial problem to identify the correct one: the tunneling method has made the choice. 
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